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The specific role played by small gold nanoparticles supported on the rutile 
TiO2(110) surface in the processes of adsorption and dissociation of H2 is dis- 
cussed. It is demonstrated that the molecular and dissociative adsorption of H2 on 
Au„ clusters containing n = 1, 2, 8 and 20 atoms depends on cluster size, geometry 
structure, cluster flexibility and the interaction with the support material. Rutile 
TiO2(110) support energetically promotes H2 dissociation on gold clusters. It is 
demonstrated that the active sites towards H2 dissociation are located at corners 
and edges on the surface of the gold nanoparticle in the vicinity of the support. The 
low coordinated oxygen atoms on the TiO2(110) surface play a crucial role for H2 
dissociation. Therefore the catalytic activity of a gold nanoparticle supported on 
the rutile TiO2(110) surface is proportional to the length of the perimeter interface 
between the nanoparticle and the support. 



1 Introduction 

Since the pioneering work of Haruta on the oxidation of carbon monoxide by 
small gold nanoparticles supported by metal oxides, ^ an extensive interest has 
been devoted to understanding the catalytic properties of gold. Such an interest is 
stipulated by the fact that gold nanoparticles are active even at room temperatures 
that makes them unique catalysts for many industrial applications.*^ 

^On leave from: V. A. Fock Institute of Physics, St Petersburg State University, 198504 St 
Petersburg, Petrodvorez, Russia. 
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The most explored type of catalytic reactions with gold nanoparticles are re- 
actions of oxidation and epoxidation, including the oxidation of carbon monoxide 
at mild temperatures, alcohol oxidation, the direct synthesis of hydrogen peroxide 
and alkene epoxidation; see, e.g., refs. 111411101 and references therein. In spite 
of intensive theoretical and experimental studies the origin of catalytic activity of 
gold remains highly debated. 

It is commonly accepted that several factors can influence the catalytic activ- 
ity of gold. The most important factor is the size effect. It has been shown that 
the unique properties of gold in oxidation reactions emerge when the size of cat- 
alytic particles decreases down to 1-5 nm or even le ss; whi le larger sized particles 
and the bulk form of gold are catalytically inactive.^MlllI n h^s been shown that 
small gold clusters consisting of a few atoms can also possess extraordinarily high 
catalytic activity. ^^^^ On the one hand, the size effects in gold nanocatalysis are 
determined by quantum effects, resulting from the spatial confinement of the va- 
lence electrons in the cluster; on the other hand, in such clusters a dominant 
fraction of atoms are under-coordinated (in comparison with the bulk), hence they 
exhibit an enhanced chemical reactivity. ^^-^^ 

Interaction with the support material is another important factor that consid- 
erably influences the chemical reactivity of the gold nanoparticles. ^^^^^^ Most 
experimental studies on the catalytic properties of gold have been performed for 
gold nanoparticles supported on the surfaces of metal oxides; see, e.g., refs. I4I19I 
and references therein. It has been demonstrated that the catalytic activity of gold 
nanoparticles depends on the presence of defects in the support material (e.g. F- 
center defects), charge transfers from the support, the chemical state of the sur- 
face, doping, moisture, special additives, and the presence of adsorbates on the 
cluster surface, including the reactant molecule itself; see, e.g., refs. 41101161191 - 
|24land references therein. However, a work by Turner et al. presents strong exper- 
imental evidence that small gold entities (~ 1.4 nm) derived from AU55 clusters 
and deposited on an inert support can also be efficient and robust catalysts for ox- 
idation reactions. Hence, the interaction with the support is important but not the 
determining factor in catalytic activity of gold; therefore even free clusters can be 
effective catalysts. 

Heterogeneously catalyzed hydrogenation is another type of reactions where 
gold nanoparticles have shown their great potential as catalysts. It has been 
demonstrated experimentally that gold nanoparticles supported on metal oxides 
such as Si02, AI2O3, Ti02, ZnO, Zr02, and Fe203 are effective catalysts for 
selective hydrogenation of several classes of organic molecules, including a, B- 
unsaturated aldehydes, unsaturated ketones, and unsaturated hydrocarbons. l^^^ 
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Moreover, supported gold nanoparticles are very selective for the direct formation 
of hydrogen peroxide from H2/O2 mixtures.^ Mechanisms of such catalytic reac- 
tions are largely not understood. Experimental results demonstrate that molecular 
hydrogen does not bind to the clean gold extended surface, and molecular ordis- 
sociative hydrogen adsorption occurs on the supported gold nanoparticles. l^^^ It 
has been suggested that hydrogen reacts and dissociates on low coordinated gold 
atoms in the corner or at edge positions on the surface of gold nanoparticles. 
It has been also shown that the shape of the gold particles plays a role on their cat- 
alytic activity. However, in a recent work by Haruta et al it has been found that 
the catalytic activity of the gold nanoparticles supported on the rutile TiO2(110) 
surface depends on the number of gold atoms located at the perimeter interface 
of the supported gold nanoparticles, irrespective of the cluster size. It was also 
supposed that active sites for H2 dissociation may be formed by a combination of 
gold atoms and oxygen atoms from Ti02 at nanoparticle-surface interface. 

Surprisingly, very little attention has been paid to theoretical investigations 
of the hydrogenation reactions on gold nanoparticles. Interaction of molecular 
hydrogen with gold atoms and small gold clusters has been theoretically studied 
in refs. 38l-l45l It has been demonstrated, that the flexibility of the cluster struc- 
ture plays a key role in the bonding and dissociation of H2 on Auj4 and Au29.'^ 
Furthermore, authors in ref. 43 have concluded that dissociation of H2 is only 
possible when a cooperation between four active low coordinated Au atoms on 
the cluster surface is allowed. However, in ref. |46| the compelling evidence has 
been presented that H2 adsorbs and dissociates with small activation barriers on 
low coordinated Au atoms situated in corner positions of different AU25 nanopar- 
ticles without further conditions. Thus, there is no need for cooperation between 
several active Au atoms to dissociate H2 and dissociation can occur spontaneously 
on a single low coordinated Au atom, when coordination number is four. 

There are relatively few theoretical studies on the catalytic hydrogenation re- 
actions on gold clusters. These include wo rks on formation of hydrogen peroxide 
from H2 and O2 over small gold clusters acrolein hydrogenation on AU20 
nanoparticle;*^ and combined experimental and theoretical investigation of the 
unusual catalytic properties of gold nanoparticles in the selective hydrogenation 
of 1,3-butadiene toward cis-2-butene.^^ 

Despite of the clear experimental evidence indicating a strong influence of the 
support materials on the processes of adsorption and dissociation of molecular 
hydrogen and related hydrogenation processes on the supported gold nanoparti- 
cles, ^^^^ there are no theoretical studies concerning the role of the support with 
an exception of recent works by Boronat et al.^^ and Florez et. al.^^ In ref. iST] 
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authors perfomied an elegant theoretical study on elucidation of active sites for 
H2 adsorption and activation on Auj3 cluster supported on the anatase Ti02 (001) 
surface. It has been shown that Au atoms that are active for H2 dissociation must 
have a net charge close to zero, be located at corner or edge low coordinated po- 
sitions, and not be directly bonded to the support. Auj3 isomers consisting 
of two layers of gold atoms have the largest number of potentially active sites 
where H2 can be adsorbed and activated. These active atoms are located on the 
top layer of the Auj3 cluster that does not have contact with the support.'^ It has 
been demonstrated that the presence of O vacancy defects on the anatase Ti02 
(001) surface can stabilize the most active two-layer Auj3 particles. 1^ In ref. l52lit 
was shown that small two-dimensional Au nanoparticles supported on TiC(OOl) 
can dissociate H2 in a more efficient way than when supported on oxides. The 
active sites for H2 dissociation on Au/TiC(001) are located at the particle edge 
and in direct contact with the underlying substrate. 

In the present paper we report the results of a systematic theoretical study of 
the adsorption and dissociation of H2 on free and TiO2(110) supported Au„ clus- 
ters containing n = 1, 2, 8 and 20 atoms. In our work we want to clarify how 
hydrogen molecule adsorbs and dissociates on the gold nanoparticles. What fac- 
tors are the most important for H2 dissociation on the surface of gold? Does the 
dimensionality of the cluster play a role in H2 adsorption and dissociation? What 
is the role of the support? Is it possible to promote dissociation of H2? Whether 
or not the charge transfer between the adsorbed H2 and the gold cluster plays any 
role in H2 catalytic activation and dissociation? Why the catalytic activity of gold 
nanoparticles supported on the rutile TiO2(110) surface depends on the number 
of gold atoms located at the perimeter interface, irrespective to the particle size 
in the case of H2 dissociation, but it depends on cluster size in the case of O2 
dissociation? Finally, what are the important directions for further theoretical in- 
vestigation of the considered processes? Answering these questions is the central 
aim of the present work. 



2 Computational Details 

The calculations are carried out using density-functional theory (DFT) with the 
gradient-corrected exchange-correlation functional of Perdew, Burke and Er nzer- 
hof (PBE).^^ The atom-centered, strictly confined, numerical basis functions 
are used to treat the valence electrons of H, Au, Ti and O atoms with the reference 
configurations ls^2p^?>(PAf for H, es^ep^Sd^'^Sf for Au, As^Ap^?>(flAf for Ti, 
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and 2s^2p^3d^4f^ for O. Double-^ plus polarization function (DZP) basis sets are 
used for Au, Ti and O, and triple-^ plus polarization function (TZP) for H. The re- 
maining core electrons are represented by the Troullier-Martins norm-conserving 
pseudopotentials^^ in the Kleinman-Bylander factorised form.^^ Relativistic ef- 
fects are taken into account for Au. 

Calculations have been carried out with the use of the SIESTA code.'^^^ Pe- 
riodic boundary conditions are used for all systems, including free molecules and 
clusters. In the latter case the size of a supercell was chosen to be large enough to 
make intermolecular interactions negligible. 

Within the SIESTA approach, the basis functions and the electron density are 
projected onto a uniform real- space grid. The mesh size of the grid is controlled 
by an energy cutoff, which defines the wavelength of the shortest plane wave 
that can be represented on the grid. In the present work the energy cutoff of 
200 Ry is chosen to guarantee convergence of the total energies and forces. The 
self-consistency of the density matrix is achieved with a tolerance of 10~^. For 
geometry optimization the conjugate-gradient approach was used with a threshold 
of 0.02 eV A-i. 

The range of the basis pseudoatomic orbitals (PAO) is limited by an energy 
shift parameter which defines the confinement of the basis functions. It has been 
shown that systema tic co nvergence to physical quantities can be obtained by vary- 
ing the energy shift.'^^^ A common energy shift of 50 meV is applied for Au, Ti 
and O on the basis of careful comparison of the obtained theoretical results with 
the available reference data for the free Au2 dimer, small gold clusters, bulk Au 
and Ti02. Thus the calculated values of the dissociation energy and bond length 
in Au2 (2.29 eV, 2.572 A) are in an excellent agreement with those of earlier 
experimental studies, (2.31 eV, 2.472 A).*^The optimized lattice constant of the 
face-centered cubic (fee) structure of bulk gold is 4.208 A, which is 3% larger than 
the experimental value of 4.079 A.^^ This is a general feature of PBE functional 
to overestimate the lattice constants for various types of solids.'^ In addition we 
tested the ability of our approach to reproduce the optimized structures and isomer 
sequences of small free gold clusters consisting of up to 10 atoms. The obtained 
structures are in a good agreement with those reported in our recent work ref. l24l 
and previous theoretical studies; see, e.g., refs. l64l - [68l 

Basis set for hydrogen was optimized with the use of the Nelder-Mead sim- 
plex method^ according to the procedure described in ref. |55] The dissociation 
energy, D^, and bond length in H2 (4.55 eV, 0.750 A) are in a good agreement 
with experimental data (4.74 eV, 0.741 A).^^ We also tested the gold - hydrogen 
interaction optimizing AuH dimer. The calculated dissociation energy and bond 
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length in AuH (3.09 eV, 1.547 A) are in a good agreement with experimental data 
(3.36 eV, 1.524 A).^^ 

The rutile structure belongs to the P42/mnm tetragonal space group. The 
primitive unit cell contains two Ti02 units. The unit cell parameters determined 
by neutron diffraction are found to be a=b=4.587 A, c=2.954 A at 15 
In the present work the rutile lattice was optimized using the Monkhorst- 
6x6x9 k-point mesh for Brillouin zone sampling. The calculated lattice param- 
eters a=4.594 A and c=2.959 A are in excellent agreement with the experimental 
values. 

The optimized lattice of the bulk rutile was used to construct slabs for the 
TiO2(110) surface. The six-layer slab containing four units of Ti02 represents the 
element of the (110) rutile face with the surface area of \/2axc. In the present 
work we study adsorption of Au„ clusters containing n = 1, 2, 8 and 20 atoms 
on the rutile (110) surface. In the case of Au and Au2 the TiO2(110) surface is 
modeled by the p(2 x 5) slab (40 units of Ti02 per slab), while for Aug and AU20 
use p(3x6) (72 units of Ti02) and p(4x9) (144 units of Ti02) slabs, respectively. 
The periodically replicated slabs are separated by the vacuum region of 25 A in 
the (110) direction. It has been shown that the surface energy, the energy gap and 
interlayer distances in thin films of rutile oscillate with decreasing amplitude as 
the number of layers increases. 1^^^^ The choice of the six-layer slab provides good 
convergence of the obtained results at moderate computational cost. Similar slab 
model has been used in ref . |76l to study an O2 supply pathways in CO oxidation 
on Au/TiO2(110). 




Figure 1: Side view of the TiO2(110) p(2x5) slab. The bridge-site twofold coor- 
dinated 0(2) atom, threefold coordinated 0(3) atom, five- and sixfold coordinated 
Ti(5) and Ti(6) atoms are marked, respectively. Relaxation of the TiO2(110) sur- 
face results in a surface rumpling. 
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Fig. [T] presents the side view of the TiO2(110) p(2x5) slab. The TiO2(110) 
surface contains atoms with different coordination. The fivefold coordinated Ti 
atoms are undercoordinated relatively to the bulk sixfold coordinated Ti. The O 
atoms raised above the plane of Ti atoms are twofold coordinated. These atoms 
are called bridge- site O. The oxygen atoms lying in the plane of Ti atoms are 
threefold coordinated. In all calculations the bottom two layers in the slabs are 
fixed, and all other atoms are fully relaxed. Relaxation of the TiO2(110) surface 
results in a surface rumpling, as it is shown in Fig. [T] The undercoordinated Ti and 
O atoms are moving inward by 0.10 A, decreasing the distance with subsurface 
atoms. The sixfold coordinated Ti atoms and threefold coordinated O atoms are 
moving outward by 0.1 A and 0.2 A, respectively. The similar puckered structure 
of the rutile (110) surface has been reported in refs. i73l75l 

3 Results 

3.1 Hydrogen adsorption on free gold clusters 

The molecular and dissociative adsorption of H2 on the free gold clusters of dif- 
ferent sizes, chains of gold atoms, and extended gold surfaces have been the sub- 
ject of a number of theoretical studies.^^^^ ^^^^"^^ It was demonstrated that 
H2 effectively adsorbs and dissociates at the low coordinated Au atoms regard- 
less if they are in g old cl usters or at extended line defects, such as monoatomic 
rows on Au surface.l^^^ However, in the case of small Au clusters, coordination 
alone cannot explain all the features of the reactivity of gold clusters for H2 dis- 
sociation.^^ Several major factors are determining the catalytic properties of free 
clusters: the structure associated with individual atoms (i.e. coordination), the 
binding structure, and the border effects.^ It was demonstrated that the planar 
gold structures are not reactive, when the H2 approaches the cluster from the top. 
However, dissociation of H2 can be promoted when planar clusters are folded or 
the H2 approaches clusters in their plane. The flexibility can play an important 
role in H2 dissociation, because of the large atomic displacements in the cluster.!^ 
Finally, the size, the structure and the charge state of gold clusters influence on 
the molecular and dissociative adsorption of H2.'^ 

In the present work we consider adsorption of H2 on gold clusters Au„ con- 
sisting of ;i = 1, 2, 8, and 20 atoms in order to elucidate the size and the geometry 
effects on hydrogen adsorption. The structural properties of free gold clusters 
have been intensively studied, see, e.g., refs. I19l68l80l for a review. DFT calcula- 
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tions demonstrate that Au„ clusters favor planar, two-dimensional (2D) structures 
up to cluster size of 11-13 atoms, while Au„ of larger sizes are three-dimensional 
(3D).'^^J^21Ilxhe high stability of planar structures for small gold clusters has been 
explained by relativistic effects. The relativistic effects make the 6s and 5d atomic 
orbitals closer in energy, resulting in a sd hybridization, which in turn stabilizes 
planar configurations. The exact value of the critical size for 2D 3D transition 
in neutral Au„ is not yet clarified. Ab initio MP2 and CCSD(T) calculations per- 
formed in ref. 84 demonstrate that Aug possesses 3D structure, while the similar 
CCSD(T) calculations performed with the larger basis set favors 2D structure for 
Aug.'^ Recent CCSD(T) calculations accounting for the basis-set superposition 
error (ESSE) corrections are in favor of 3D structure for Aujq^^ in accord with 
DFT predictions. Competition of 2D and 3D structures of the small gold clusters 
is a very interesting feature that can play an important role in cluster reactivity. 
The 2D structures contain a large number of low coordinated atoms, hence they 
might possess higher catalytic activity in comparison with the 3D clusters. More- 
over, it has been shown that if the shape of neutral gold clusters is modified, it is 
possible to change their electron donor- acceptor capacities - the 2D gold clusters 
are better electron acceptors than 3D ones.^^ Since the electron transfer between 
the cluster and the adsorbate is an important driving force for the reactivity, the 
shape of the cluster can directly affect cluster's catalytic properties.!^ 

In our previous works we have optimized the most stable and isomer geome- 
tries of small gold clusters consisting up to 10 atoms within DFT B3PW91/LANL2DZ 
approach. In the present work we have re-optimized the earlier obtained ge- 
ometries of Au„ using the PBE/DZP SIESTA method. The obtained structures are 
in an excellent agreement with those reported in previous theoretical studies; see, 
e.g., refs. I24l64» 68l 

In Fig. |2| we present optimized geometries of Au„ systems in the case of 
molecular and dissociative adsorption of H2 on free Au„ clusters. In order to ob- 
tain the most stable configuration of H2 adsorbed on Au„, we have created a large 
number of starting geometries by adding H2 molecule in different positions (up 
to 30) on the surface of the considered clusters. The starting structures have been 
optimized further without any geometry constraints. We have successfully used 
a similar approach to find the optimized geometries of O2 and C2H4 molecules 
adsorbed on small neutral, anionic and cationic gold clusters. ^^^^^^^ In addition 
to the most stable 2D structure of Aug and 3D structure of AU20 consider the 
closest in energy 3D isomer of Aug and 2D isomer of AU20 order to under- 
stand how dimensionality and structural properties affect the cluster reactivity. 
The calculated relative energies for the higher isomeric Aug (3D) and AU20 (^^) 
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Au 

Au2 



AUb (2D) 



AUb (3D) 
(+0.18 eV) 

AU20 (3D) 



Au£o (2D) 
(+1.95 eV) 




Figure 2: Optimized geometries of Au„ clusters in the case of molecular 
(left) and dissociative (right) adsorption of The numbers in the parentheses 
correspond to the relative energy for the higher isomeric free Aug (3D) and AU20 
(2D) clusters with respect to the corresponding Aug (2D) and AU20 (^^) lowest 
energy structures. 
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structures with respect to the corresponding lowest energy structures are 0.18 eV 
and 1.95 eV, respectively. The binding energies for molecular, E^^^ , and dissocia- 
tive, Ep\ adsorption of H2 on free Au„ as well as the H— H bond length r^l^i in 
H2 adsorbed molecularly are summarized in Table [T] 

Table 1: Binding energies for molecular, E^^^ , and dissociative, Ep\ adsorp- 
tion of H2 on free Au„ clusters and the H— H bond length, r^^^' ^^^^ 
molecular adsorption of 



Cluster 


rg-'n (A) 




E'l- (eV) 


Au 


0.780 


0.13 


0.04 


AU2 


0.845 


0.59 


0.59 


Aug (2D) 


0.806 


0.26 


0.68 


Aug (3D) 


0.788 


0.13 


0.52 


AU20 (3D) 


0.768 


0.09 


0.14 


AU20 (2D) 


0.762 


0.09 


0.72 



The binding energy of H2 adsorbed on a free Au„ cluster is defined as 

^^oi4is ^ ^ ^^^^(j^^) _ £-W«(H2_AuJ, (1) 

where £'^'^J^^'^'^(H2— Au„) denotes the total energy of the compound system, H2— Au„, 
while Etot{^^n) Etot{)^2) ^^^^^ energies of the non-interacting frag- 

ments Au„ and H2, respectively. 

Results of our calculations demonstrate, that H2 binds weakly to Au as a 
molecule, with £"^^^=0.13 eV. Dissociation of H2 on Au is not favorable ener- 
getically. The H— H bond length in H2 adsorbed on Au is slightly enlarged in 
comparison with the free H2. Adsorption of H2 on AU2 is relatively strong, with 
the binding energy equal to 0.59 eV both for the molecular and the dissociative 
adsorption. Table [I] demonstrates that the H— H bond length increases to 0.845 A 
in H2— Au2, indicating that H2 is highly activated. 

With the further increase in cluster size the binding energy calculated for 
molecular adsorption of H2 decreases to 0.26 eV (0.1 3eV) for Aug (2D) (Aug 
(3D)) and 0.09 eV both for 3D and 2D isomers of AU20. other hand, the 

dissociative adsorption of H2 becomes energetically favorable both for Aug and 
AU20. binding energy of H2 on free gold clusters strongly depends on the 
cluster's isomer state. Thus, dissociative adsorption of H2 on 2D structures of 
Aug and AU20 energetically favorable compared with the corresponding 3D 
structures. Therefore, dissociation of H2 on small gold clusters can be promoted 
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if we can find a way to enhance the stabiUty of 2D structures for n > 13. That 
should be possible for supported clusters, where an attractive interaction with the 
surface stabihzes oblate configurations.^^ Thus, varying the interaction with the 
support one can change the cluster shape, and hence it can be possible to tune its 
reactivity. Another factors that influence H2 dissociation are coordination of Au 
atoms interacting with hydrogen and flexibihty of cluster structure. It is seen from 
Fig. |2|that H2 dissociates at the low coordinated comer Au atom with formation 
of the slightly bended H— Au— H bond. It is worth to note, that the Au— Au bond 
length in AU2 does not change noticeably upon molecular adsorption of but 
increases up to 2.753 A for dissociative adsorption. In the case of Aug and AU20 
dissociation of H2 on the cluster surface is accompanied by the structural (at least 
local) transformations. Therefore we can confirm conclusion made in ref. 43 that 
the structural flexibility of the gold clusters might be an important factor for H2 
dissociation. 

3.2 Hydrogen adsorption on Au^/TiO2(110) 

How does the interaction with the support influence the reactivity of gold clus- 
ters? Can a supported nanoparticle exhibit novel properties that are not found 
for free particles? Recent experiments performed by Haruta et al clearly demon- 
strate that the perimeter interface of the gold nanoparticles supported on the rutile 
Ti02(l 10) surface plays a crucial role in the process of hydrogen dissociation.^-^ It 
was suggested that the active sites for H2 dissociation might be formed by a com- 
bination of gold atoms and oxygen atoms from Ti02 at the nanoparticle-surface 
interface.'^ 

In the present work we perform a systematic theoretical study of the structure 
and energetics of H2 adsorbed on Au^, Au2, Aug (2D), Aug (3D), AU20 (^^) 
AU20 (^^) clusters supported on the rutile TiO2(110) surface. The analysis of the 
energetics of H2 adsorption on the supported gold clusters provides deep insights 
into the nature of bonding and dissociation of the H2 molecule and reveals the 
details of the reaction mechanism. 

As discussed above the dissociation of H2 on free gold clusters is governed by 
the interplay of structural factors, coordination of the adsorption center and cluster 
flexibility. The complexity of the changes in structural and electronic features of 
gold nanoparticles becomes more diverse when one deposit nanoparticle on the 
support. 

In order to obtain the most stable geometries of Au„ supported on the rutile 
TiO2(110) surface we have created a large number of starting configurations by 
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adding Au„ isomers with different orientation in respect to the surface. These 
structures have been fully optimized on the rutile surface, with accounting for the 
relaxation of all the gold atoms as well as the top four layers of the six layer slab 
of the rutile. The bottom two layers in the slab were fixed. Thus, we have taken 
into account deformations of the cluster structure due to the interaction with the 
support as well as structural relaxations on the rutile TiO2(110) surface due to its 
interaction with the supported cluster. 




Eb= 0.69 eV Eb= 1 .65 eV E^= 1 .47 eV 




Eb= 1 .32 eV Eb= 1 .64 eV E^= 2.53 eV 



Figure 3: The most stable geometries of Au„ clusters (n = I, 2, S, and 20) opti- 
mized on the rutile Ti02(l 10) surface. Top views of Au and Au2 on Ti02(l 10) are 
shown in the inserts. The binding energies of Au„ on Ti02(l 10) are marked below 
the corresponding structures. Numbers in parentheses correspond to the relative 
energy for the higher isomeric free Aug (3D) and AU20 (^^) clusters with respect 
to the corresponding free Aug (2D) and AU20 (•^^) lowest energy structures. 

Fig. [3]presents the most stable geometries of Au„ clusters (n= 1, 2, 8, and 20) 
optimized on the rutile Ti02(l 10) surface. We found that gold tends to maximize 
interaction with the 0(3) and Ti(5) atoms on TiO2(110) surface. Thus, Au atom 
tends to occupy a position above the row of 0(3) surface atoms, simultaneously 
bridging two 0(3) atoms in the direction parallel to the 0(3) row, as well as Ti(5) 
and 0(2) atoms, in the direction perpendicular to the 0(3) row. In the case of Au2, 
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both of Au atoms occupy positions above of the rows of 0(3) atoms, symmetri- 
cally in respect to the row of Ti(5) atoms. The Au— Au bond length of 2.57 A in 
Au2 matches the distance between the two rows of 0(3) atoms (2.53 A ). The ge- 
ometry structures of two- and three dimensional isomers of Aug on the Ti02(l 10) 
surface are slightly deformed if compared with the free gold clusters. The Aug 
(2D) cluster maximizes the interaction between the cluster edge and the row of 
0(3) atoms on the TiO2(110) surface. The angle between the cluster's plane and 
the rutile surface is 26°. The supported Aug (3D) isomer prefers to maximize 
its interaction with the rutile TiO2(110) surface, orienting the 5 -atom face paral- 
lel to the support, as shown in Fig. [3] The binding energies of Aug (2D) and 
Aug (3D) isomers on TiO2(110) are 1.65 eV and 1.64 eV, respectively. A tetrahe- 
dral Au2o (^^) cluster binds to the surface, maximizing interaction between the 
cluster edge and rows of 0(3) and Ti(5) surface atoms; while AU20 (^^) isomer 
binds by two its edges to two parallel rows of Ti(5) surface atoms, slightly bend- 
ing its planar structure, to avoid interaction with the row of 0(2) bridge atoms in 
the middle. The planar AU20 structure binds to the surface much more strongly 
than the 3D structure, because it has a wider contact area. Although AU20 (•^^) 
structure is still energetically favorable when adsorbed on the TiO2(110) surface 
in comparison with the 2D structure, the energy difference between two configu- 
rations decreases from 1.95 eV for free clusters to 0.89 eV for supported clusters, 
respectively. If we could further increase the interaction with the support, we 
would make 2D structure energetically favorable. Such an effect has been re- 
ported for AU20 clusters deposited on thin MgO(lOO) films supported on Mo(lOO) 
substrate. By changing the thickness of the metal-oxide film, one can tune the 
interaction of the cluster with the support, and thus increase its wetting propensity 
and induce a dimensionality crossover from 3D cluster structures on MgO(lOO) 
to the energetically favored 2D geometries on the metal- supported films .'^ 




Figure 4: Optimized geometries for molecular (left) and dissociative (right) ad- 
sorption of H2 on the pure rutile Ti02(l 10) surface. 
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Fig. |4| presents the most stable geometries calculated for the molecular and 
dissociative adsorption of H2 on the pure rutile TiO2(110) surface. Our calcula- 
tions demonstrate that the hydrogen molecule adsorbs on the pure rutile Ti02(l 10) 
surface on top of the Ti(5) atom with the binding energy of 0.14 eV. The H— H 
bond length of 0.78 A is slightly increased in comparison with the free H2. The 
dissociative state of H2 on TiO2(110) corresponds to the situation when both H 
atoms form the OH group with the low coordinated 0(2) bridge atoms on the ru- 
tile surface. The binding energy of the dissociated configuration of H2 is 1.50 
eV; however, the distance between two rows of 0(2) atoms on TiO2(110) is 6.50 
A, which is too large to promote dissociation of the H2 adsorbed on top of Ti(5) 
atom. In this case H2 would dissociate in the vicinity of the adsorption center as a 
first step, followed by adsorption of H atoms on 0(2). This process would require 
to overcome a dissociation barrier much higher than the energy of the molecular 
adsorption. In this situation H2 would likely desorb from the surface and fly away 
rather than dissociate. Therefore, to promote H2 dissociation on TiO2(110), it is 
necessary to have adsorption centers on the rutile surface in the vicinity of low 
coordinated 0(2) atoms. Supported gold clusters can serve as a source of such 
centers. 




Figure 5: Optimized geometries for molecular (left) and dissociative (right) ad- 
sorption of H2 on Au/Ti02. 

Fig. |5]presents optimized geometries calculated for molecular and dissociative 
adsorption of H2 on Au/Ti02. It is seen from Table |2| that H2 adsorbs molecularly 
on Au/Ti02 with a binding energy of 1.15 eV, which is considerably larger than 
the corresponding energy calculated for H2 adsorbed on free Au atom. The ad- 
sorbed H2 is highly activated, with the H-H bond length r^^^p^=0.906 A. As a 
result of H2 adsorption, the supported Au atom shifts slightly towards the row 
of the low coordinated 0(2) bridge atoms. In the previous section we saw that 
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Table 2: Binding energy calculated for molecular, E^^, and dissociative, E^^\ 
adsorption of H2 on Au„/Ti02, and the H— H bond length, v^Iy^, in the case of 
molecular adsorption of 



Cluster 


rg^^H (A) 




Ei- (eV) 


Au 


0.906 


1.15 


2.06 




0.821 


0.37 


0.73^^) 0.90^^^) 1.55^^) 


Aug (2D) 


0.792 


0.24 


0.75^^) 1.55(^) 


Aug (3D) 


0.753 


0.09 


0.55^^) 0.95(^) 


AU20 (3D) 


0.795(^) 0.806^^') 


0.17(^) 0.13^^') 


0.20^^) 0.83(^) 


AU20 (2D) 


0.798^^) 0.805^^') 


o.ii(^) o.ii(^') 


0.55(^) 1.42(^) 



H2 dissociation is not favorable on the free Au atom. However, in the case of 
the supported Au atom, dissociation of H2 can occur with formation of the OH 
group with the 0(2) atom located either in the nearest to Au row of 0(2) or in the 
next row of 0(2), as is shown in Fig. |5] In the latter case the calculated binding 
energy E^^^=2.06 eV is higher, not only if compared with the binding energy for 
molecular adsorption of H2 on Au/Ti02(l 10), but also if compared with the bind- 
ing energy of the dissociated state of H2 on the pure TiO2(110) surface. Thus, 
we can conclude, from the energetic point of view, that the TiO2(110) support 
considerably promotes H2 dissociation on Au atom. 

The optimized geometries calculated for molecular and dissociative adsorption 
of H2 on Au2/Ti02 are shown in Fig. [6| Figs. Qa) andQb) demonstrate that 
H2 binds to the supported Au2 in a similar way as for the free Au2. However 
the binding energy calculated for the molecular adsorption of H2 on Au2/Ti02 is 
lower than the corresponding energy obtained for free AU2, while the dissociated 
configuration of H2 on Au2/Ti02 is more stable if compared with dissociative 
adsorption of H2 on the free AU2. Therefore, the interaction of AU2 with the 
support results in the energetic promotion of H2 dissociation on the supported 
Au2. We found, however, that the geometry configuration shown in Fig. Qb) is 
not the most stable one for H2 dissociation. We found that the hydrogen atom can 
migrate to the row of 0(3) atoms, or to the row of 0(2) bridge atoms on the rutile 
surface, forming OH bonds, as it is shown in Figs. [6fc) and[6fd), respectively. 
Hydrogen atom always binds more strongly with the low coordinated 0(2) atoms, 
as can be seen from Table |2l 

Figs. |7]^a),|7][b),|7]^d) and|7][e) demonstrate that optimized geometries of H2 in 
the case of its molecular and dissociative adsorption on the supported Aug (2D) 
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(a) (b) 




(c) (d) 



Figure 6: (a) Optimized geometry in the case of molecular adsorption of H2 on 
Au2/Ti02. Optimized geometries in the case of dissociative adsorption of H2 on 
Au2/Ti02: (b) dissociation of H2 on supported AU2; (c) dissociation of H2 with 
formation of the OH group on 0(3) surface atom; (d) dissociation of H2 with 
formation of the OH group on 0(2) bridge atom. 
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Figure 7: Optimized geometries in the case of (a) molecular adsorption of H2 
on Aug(2D)/Ti02; dissociative adsorption of H2 on Au2(2D)/Ti02: (b) dissocia- 
tion of H2 on supported Au2(2D); (c) dissociation of H2 with formation of the OH 
group on 0(2) bridge atom. Optimized geometries in the case of (d) molecular ad- 
sorption of H2 on Au3(3D)/Ti02; dissociative adsorption of H2 on Au2(3D)/Ti02: 
(e) dissociation of H2 on supported Au2(3D); (f) dissociation of H2 with formation 
of the OH group on 0(2) bridge atom. 
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and Aug (3D) clusters similar to those obtained for the corresponding free Aug 
clusters. Note that the dissociative adsorption of H2 on Aug (3D) results in the 
considerable rearrangement of the gold structure as shown in Fig. |7f e). Tables [T] 
and |2| demonstrate that the binding energies calculated for molecular adsorption 
of H2 on free and supported Aug (2D) and Aug (3D) are quite similar. However, 
the dissociative configuration of H2 on the supported Aug (2D) and Aug (3D) 
clusters is slightly promoted energetically. Migration of H on the low coordinated 
0(2) atom as it is shown in Figs. |7fc) and |7ff), results in formation of the OH 
bond and considerable increase in binding energies calculated for H2 adsorbed 
dissociatively. It is seen from Fig. |7]^c) that the planar structure of Aug undergoes 
drastic rearrangements upon H2 dissociation. The dissociative adsorption of H2 on 
Aug is energetically favorable if compared with the adsorption on the 3D isomer 
of Aug. 

The AU20 (^^) ^^20 (^^) clusters supported on TiO2(110) are the largest 
systems studied in the present work. In the previous section it was shown that the 
hydrogen molecule binds weakly to the free AU20 (•^^) ^^20 (^^) clusters. 
We have found that, in the case of the supported AU20, hydrogen molecule 
adsorbs on top of the Ti(5) atoms on the rutile surface in the vicinity of AU20, 
rather than directly on the AU20 clusters, as it is shown in Figs. [8]^a),|8]^d) and[8]^d'). 
The geometry structure when H2 adsorbs at the vertex of Au2o(3D), bridging Au 
vertex atom in the cluster and the low coordinated 0(2) atom on the rutile surface 
is also stable (Fig. [8]^a')). Further H2 dissociation can occur either at the vertex 
of AU20 (•^^) acute-angled corner of AU20 (2^)' shown in Fig. |8][b) 

and Fig. [8]^e), respectively. H2 dissociation at the acute-angled corner of AU20 
(2D) is accompanied by severe local structural rearrangement of Au atoms at the 
corner of AU20 ('^^) migration of one H atom to the 0(2) atom on the surface. 
Note that the dissociation of H2 at the obtuse-angled comer of the supported AU20 
(2D) cluster is not energetically favorable, in comparison to the case of free AU20 
(2D). The most stable configurations of the dissociated H2 are those where one 
of the H atoms binds to the edges of AU20 (•^^) ^^20 ('^^) ^^^^ oriented 
perpendicularly to the rows of 0(2) surface atoms, while another H forms the OH 
bond with the 0(2) atom on the rutile surface, as is shown in Figs. [8fc) and[8ff). 
We found several geometrical configurations of this type with binding energies of 
1.19 - 1.42 eV. Dissociation at the edge of the supported 2D isomer of AU20 
always energetically favorable. 

Thus, combination (interplay) of several factors such as geometry structure, 
cluster dimensionality, presence of the low coordinated oxygen atoms in the vicin- 
ity of the cluster- surface interface, etc. are important for H2 dissociation. 
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(c) (f) 



Figure 8: Optimized geometries for H2 adsorbed on Au2o(3D)/Ti02 (left) and 
Au2o(2D)/Ti02 (right), (a) and (a') Molecular adsorption of H2 in the vicinity 
of Au2o(3D)/Ti02; (b) dissociative adsorption of H2 at the vertex of Au2o(3D) 
with formation of the OH group; (c) dissociative adsorption of H2 at edge of 
Au2o(3D) with formation of the OH group; (d) and (d') molecular adsorption of 
H2 in the vicinity of Au2o(2D)/Ti02; (e) dissociative adsorption of H2 at the corner 
of Au2o(2D) with formation of the OH group; (f) dissociative adsorption of H2 at 
edge of Au2o(2D) with formation of the OH group. 
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It is well known that, in the case of catalytic oxidation reactions by molecular 
oxygen on free and supported gold nanoparticles, the charge transfer from the gold 
to the antibonding orbital of O2 is responsible for the catalytic activation and dis- 
sociati on of Q9. However, in the case of H2 dissociation the analysis of the Bader 
chargesl^^^ demonstrate that there is no considerable charge transfer between the 
adsorbed hydrogen and the gold clusters (free or supported). We found the largest 
charge transfer occurs for H2 adsorbed molecularly on Au/Ti02 (Fig. |5]). In this 
case the calculated Bader charge localized on H2 is +0.11e, where e is the ele- 
mentary charge. Although such a charge transfer is relatively small, it might be 
responsible for the strong enlargement of the H— H bond length in H2— Au/Ti02 
up to 0.906 A. Nevertheless, we have not found the direct correlation of the hy- 
drogen dissociation with the charge transfer between hydrogen and gold atoms. 
On the other hand, formation of the OH group with the bridge 0(2) atom on the 
rutile surface is accompanied by the large charge transfer from H atom to 0(2), 
resulting in a Bader net charge of H in OH equal to -i-0.75e. Thus, we can conclude 
that the catalytic activity of gold nanoparticles for O2 dissociation would depend 
on the electronic structure and the size of the nanoparticles; however in the case of 
H2 dissociation it will be proportional to the number of gold atoms located in the 
vicinity of the low coordinated 0(2) atoms at the nanoparticle - surface interface. 



4 Conclusions 

The present theoretical study demonstrates that adsorption of H2 on Au„ (n= 1, 2, 
8, 20) strongly depends on cluster size, geometry structure, flexibility and inter- 
action with the support material. Strong interaction with the support can stabilize 
2D structures of gold clusters on the surface. This is important finding, because 
the energy release due to H2 dissociation is the largest for the 2D isomers of 
gold clusters. Rutile TiO2(110) support energetically promotes H2 dissociation 
on gold clusters. The formation of the OH group near the supported gold clus- 
ter is an important condition for H2 dissociation. We have shown that the active 
sites towards H2 dissociation on the supported Au„ are located at comers and 
edges of the gold cluster in the vicinity of the low coordinated oxygen atoms on 
TiO2(110). Therefore catalytic activity of a gold nanoparticle supported on the 
rutile TiO2(110) surface is proportional to the length of the perimeter interface 
between the nanoparticle and the support, in accordance with the recent experi- 
mental findings. 

In the present work we have systematically studied adsorption and dissocia- 
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tion of H2 based on the assumption of the total energy criterion. Thus, we have 
identified a number of the most stable configurations for molecular and dissocia- 
tive states of H2 adsorbed on free and supported gold nanoparticles. Conditions 
leading to the stabilization of the dissociated state of H2 have been disclosed. 

However, many interesting and important questions regarding the mechanism 
of H2 dissociation on Ti02(l 10) supported Au„ clusters remain unanswered. Thus, 
it is necessary to elucidate the full reaction pathways and to estimate the reaction 
barriers for H2 dissociation. This will allow favorable reaction channels and reac- 
tion dynamics to be identified. As a next step one should consider simple hydro- 
genation reactions, such as, for example, hydrogen peroxide formation from H2 
and O2 catalyzed by the supported gold clusters. We can suggest that structures 
active for H2 dissociation, but not active enough for O2 dissociation could work 
as catalysts for H2O2 synthesis. 

In many cases we have found an ensemble of initial and final configurations 
of the adsorbed hydrogen with small differences in binding energies, but possible 
difference in the reaction barriers. Accounting for the geometry of the initial 
and final states can lead to further non-equivalence of different reaction pathways 
and may have an important impact on the change of system entropy. It has been 
arguedl^^^ that accounting for an entropy change contribution to the free energy 
of the system is necessary for correct description of the temperature dependencies 
of the branching ratios between different reaction channels, while purely energetic 
considerations based on reaction barrier heights can fail. 

Finally, we have considered the pure, defect free rutile TiO2(110) surface. 
It is well known that the surface defects, especially O vacancies in metal-oxide 
supports, or extra O atoms adsorbed on the surface can change drastically the 
catalytic activity of the supported gold clusters in oxidation reactions by molec- 
ular oxygen. ^^^^ Interaction of small gold clusters with a partially reduced rutile 
TiO2(110) support can be rather complicated; supported gold clusters can either 
receive or donate electrons to the substrate depending on the cluster size. The 
role of the surface defects in the process of H2 dissociation is not so obvious and 
requires additional systematic investigation. Such defects, in particular additional 
O atoms adsorbed on the surface in the vicinity of the cluster - surface interface 
can probably play a role in H2 dissociation and hence can modify the cluster reac- 
tivity. Understanding how to control such chemical reactions on a cluster surface 
is a vital task for nanocatalysis. 

We hope that the present discussion will stimulate further experimental and 
theoretical investigations of the processes considered. 
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